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Collinear high-gradient 0{GVfm) beam-driven wakefield methods for charged-particle acceler¬ 
ation could be critical to the realization of compact, cost-efficient, accelerators, e.g., in support 
of TeV-scale lepton colliders or multiple-user free-electron laser facilities. To make these options 
viable, the high accelerating fields need to be complemented with large transformer ratios > 2, a 
parameter characterizing the efficiency of the energy transfer between a wakefield-exciting “drive” 
bunch to an accelerated “witness” bunch. While several potential current distributions have been 
discussed, their practical realization appears challenging due to their often discontinuous nature. 

In this paper we propose several alternative current profiles which are smooth which also lead to 
enhanced transformer ratios. We especially explore a laser-shaping method capable of generating 
one the suggested distributions directly out of a photoinjector and discuss a linac concept that could 
possible drive a dielectric accelerator. 

PACS numbers: 29.20.Ej, 29.27.-a, 41.85.-p, 41.75.Fr 


I. INTRODUCTION 

In beam-driven techniques, a high-charge “drive” 
bunch passes through a high-impedance medium and ex¬ 
periences a decelerating field [1-3]. The resulting energy 
loss can be transferred to a properly delayed “witness” 
bunch trailing the drive bunch. A critical parameter as¬ 
sociated to this class of acceleration method is the trans¬ 
former ratio 


where is the maximum accelerating field behind the 
drive bunch, and E- is the maximum decelerating field 
within the drive bunch. 

Generally the transformer ratio is limited to values 
IZ < 2 due to the fundamental beam-loading theo¬ 
rem [4]. However larger values can be produced using 
drive bunches with tailored (asymmetric) current pro¬ 
files. Furthermore, it can be shown that both IZ and Ej^ 
for a given charge are maximized when the decelerating 
field over the drive bunch is constant [5]. Additionally, 
bunch current profiles that minimize the accumulated en¬ 
ergy spread within the drive bunch are desirable as they 
enable transport of the drive bunch over longer distances. 

To date, several current profiles capable of generating 
transformer ratios IZ > 2 have been proposed [5-7]. 
These include linearly ramped profiles combined with a 
door-step or exponential initial distribution [8]. More 
recently a piecewise “double-triangle” current profile 
was suggested as an alternative with the advantage of 
being experimentally realizable [9]. A main limitation 
common to all these shapes resides in their discontinuous 
character which make their experimental realization 
either challenging or relying on complicated beam- 


manipulation techniques [10, 11]. In addition these 
shapes are often foreseen to be formed in combination 
with an interceptive mask [12, 13] which add further 
challenges when combined with high-repetition-rate 
linacs [14]. 

In this paper we introduce several smooth current pro¬ 
files which support large transformer ratios and lead to 
quasi-constant decelerating fields across the drive bunch. 
We describe a simple scheme for realizing one of these 
shapes in a photoemission radiofrequecy (RF) electron 
source employing a shaped photocathode-laser pulse. Fi¬ 
nally, we discuss a possible injector configuration that 
could form drive bunches consistent with the multi-user 
free-electron laser (FEE) studied in Ref. [14]. 

II. SMOOTH SHAPES 

For simplicity we consider a wakefield-assisting 
medium (e.g. a plasma or a dielectric-lined waveguide) 
that supports an axial wakefield described by the Green’s 
function [15] 

G{C) =2Kcos{kC), (2) 

where k is the loss factor and k = 27r/A with A being 
the wavelength of the considered mode. Here C > 0 (in 
our convention) is the distance behind the source parti¬ 
cle responsible for the wakefield. In this Section we do 
not specialize to any wakefield mechanism and recognize 
that, depending on the assisting medium used to excite 
the wakefield, many modes might be excited so that the 
Green’s function would consequently consist of a summa¬ 
tion over these modes. 

Given the Green’s function, the voltage along and be¬ 
hind a bunch with axial charge distribution S{z) can be 
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obtained from the convolution [15] 


A. Linear ramp with sinusoidal head 


V{z)= r Giz-OSiOdC- (3) 

J — oo 


We take S{z) to be non vanishing on two intervals 
[0, and [^, Z] and zero elsewhere. In our convention the 
bunch head starts at 2 ; = 0 and the tail population lies 
at z > 0. We also constrain our search to functions such 
that S{z) and S'{z) = dS/dz are continuous at ^ 
Introducing the function f{z) (to be specified later), we 
write the charge distribution as 



/(^) 

noz-rm+fio 


0 


if 0 < ^ < ^, 
if ^ < z < Z, 
elsewhere. 


( 4 ) 


Based on our previous work [16] we first consider the 
following function 

f{z) = az-\-bsm{qkz), (5) 

where a and b are positive constants, k is again the spa¬ 
tial frequency seen above, and g' > 0 is an integer. Con¬ 
sequently, using Eq. 4, the axial bunch profile is written 
as 


az + bsm{qkz) 

ax H- bqk{x — <f) cos{q^k) 

-\-bsm{q^k) 

0 

V 


if 0 < ^ < ^, 

if ^ < z < Z, 
elsewhere. 


(6) 


In this section we report only on solutions pertaining to 
^ = A/2. Additional, albeit more complicated, solutions 
also exist for larger however, these solutions lead 
to additional oscillations which ultimately lowers the 
transformer ratio. 


From Eq. 3, the decelerating field then takes the form 


E-{z) = n < 


^ (oAsin^ {!f) + 

A((g^-l)(aA+27r6(-l)^g)+cos(^)(27r6g((-l)'^g^ + l)-a(g^-l)A)) 

27r2(g2 — 1 ) 


z < A/2 
z > A/2 


( 7 ) 


The oscillatory part in the tail (A/2 < z) can be elim- which leads to the following decelerating and acceler- 
inated under the condition ating fields respectively 


a (g^ — 1 ) A 

27rq ((-l)^g2 + 1) ’ 


pZ=NX 

E_^{z) = / s{z')uj{z — z')dz' 

aX^ (tt (( — 1)^ — l)g^ — 2N -1- l) + 2N) sin (27r (jV-f)) + ((-l)q2g2-l)+l) COS (27r (a — f fl) ^ 

27r2 ((-1)V + 1) • 


(8) 


E (z)- K I 2nm-l)^q^ + l) ^ ^ 

27r2((-l)‘?g2 + l) ^ 


Finally, the transformer ratio can be calculated by tak- pendix A) over the maximum decelerating field which 
ing the ratio of the maximum accelerating field (see Ap- yields 


(-i)nv-i) + i 


(11) 


Two sets of solutions occur for even and odd q which can be interpreted as a phase shift in the oscillatory part. 
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FIG. 1. Example of current profiles described by Eq. 6 
(shaded line) with the corresponding induced voltages. The 
parameters are n = 0, iV = 5, and plots (a) and (b) respec¬ 
tively correspond to the cases q — 2 and q = 3. The head of 
the bunch is at = 0. 


Additionally, larger multiples of even and odd q lead to 
more oscillations in the head which ultimately reduce the 
transformer ratio. In Fig. 1 we illustrate the simplest 
even (a) and odd (b) solutions corresponding to q = 2 
and q = 3 respectively. 


B. Linear ramp with parabolic head 

We now consider an even simpler quadratic shape 
which was inspired by our previous work [16, 17] 



FIG. 2. Example of ideal “quadratic” current profiles given 
by Eq. 13 (shaded line) with corresponding induced voltage. 
The parameters are v — and N — 3. The head of the bunch 
is at koz = 0. 


The accelerating field trailing the bunch {z > Z = NXq) 
is 

E^{z) = - [7r(2A^ - u) sm{kz) (14) 

— cos{kz)], 
yielding the transformer ratio 

7^= [l + 7r^(2A^-z/)Y/^- (15) 

In Fig. 2 we illustrate an example of the quadratic shape 
(green trace) as well as its corresponding longitudinal 
electric field (blue trace) for = 1 and = 5. 


f{z) = az^, 


which leads to the current profile 


S{z) = Z 


az 


2a^z — 
0 


if 0 < 2 ; < (f, 

if C ^ ^ ^ 

elsewhere. 


( 12 ) 


(13) 


The resulting decelerating field within the bunch is 


E_{z) = 2 k < 


r)^ sin(kz) — kz 
r)^ sin[/c(2;—^)]—sin(/c2;)-|-2/c^ 


0 


if 0 < ^ < ^, 
if C ^ ^ ^ ^5 

elsewhere. 


Again, the decelerating field can be made constant for 
2 ) G [^, Z] when ^ = i/\ with u e N. In such a case the 
previous equation simplifies to 


E_{z) = 2t<i < 


r)^ s\nikz) — kz 


k^ 


0 

\ 


if 0 < 2) < z^A, 
if z^A < 2) < Z, 
elsewhere. 


C. Comparison with other shapes 

We now turn to compare the smooth longitudinal 
shapes from the previous Section with the doorstep [8] 
and double-triangle [9] which also provide constant de¬ 
celerating fields over the bunch-tail (see Appendix B for 
our formulation of these distributions). For a fair com¬ 
parison, we stress the importance of comparing the vari¬ 
ous current profiles with equal charge. Consequently, we 
normalize each of the current profile to the same bunch 
charge 


nZ = NX 

Q= dzS{z^a)] (16) 

Jo 

where a is the scaling parameter associated with each 
bunch shape (see Section II and Appendix B), and NX 
is the total bunch length which is assumed to be larger 
than the given shape’s bunch-head length {NX > ^). For 
each distribution, the charge normalization generates a 
relationship between a and NX which enables us to rex- 
press IZ in terms of Q and a. In Tab. I we tabulate the 
analytical results for 71{N) (the conventional notation 
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[8, 9]) and lZ{Q^a)^ and also list the maximum decel¬ 
erating field for each distribution. Additionally in 
Fig. 3 we illustrate these results in a log-log plot where, 
for each distribution, the scaling parameter (a) was var¬ 
ied for a fixed charge and wavelength. To complete our 
comparison we also added the linear-ramp and Gaussian 
distributions. 

The results indicate that all of the distributions with 
constant decelerating fields over the bunch-tail dive’ on 
the same curve; additionally, by varying the scaling pa¬ 
rameter a for a distribution, you can shift a distribu¬ 
tion to have a larger (resp. smaller) IZ (resp. and 

vice-versa. Ultimately, this suggests that the distribu¬ 
tion which is simplest to make is as useful as any other 
and it can be scaled accordingly (7^, for a specific 
application. These results confirm our previous studies 
regarding the numerical investigation of the trade-off be¬ 
tween IZ and E-^ [18]. 



FIG. 3. Tradeoff curves between 7Z and EH^ for the current 
profiles listed in Tab. I. The “quadratic” and “sin ramps” 
respectively correspond to the distributions proposed in Sec¬ 
tions IIB and II A. The Gaussian and ramp distributions are 
displayed for comparison. 


III. PHOTOEMISSION OF OPTIMAL SHAPES 
VIA LASER-SHAPING 

In this section we investigate the realization of the 
quadratic distribution discussed in Section II by longi¬ 
tudinally tailoring a laser pulse impinging on a photo¬ 
cathode in a photoinjector. The resulting electron dis¬ 
tribution is then accelerated in an RF-gun and expands 
via space charge forces. If the charge density of the em¬ 
anating electron bunch is sufficiently low, the resulting 
distribution will be relativistically preserved through a 
drift; however for larger charge densities, the original lon¬ 
gitudinal distribution will morph according to the inte¬ 
grated space charge forces inside the bunch. The setup 


we consider throughout this section is depicted in Fig. 4 
and consists of a typical I + ^-cell BNL/SLAC/UGLA 
S-band RF-gun operating at 2.856 GHz surrounded by 
a solenoidal lens [19]. The large 140 MV/m) accel¬ 
eration gradients in the gun help preserve larger charge 
densities compared with e.g. L-band guns. The sim¬ 
ulations are carried with Astra [20], a particle-in¬ 
cell beam-dynamics program that includes a quasi-static 
cylindrically-symmetric space charge algorithm. The 
simulation also includes the image-charge effect which 
arises during the photoemission process, in our simula¬ 
tions the electron bunch is represented by 200,000 macro¬ 
particles. 



FIG. 4. Configuration used for the pulse-shaping simulations 
using a S-band RF gun (a). A temporally shaped laser pulse 
(b) is optimized to result in a photo-emitted electron-beam 
with current profile (c) having features similar to the distri¬ 
bution discussed in Sec. IIB. In insets (a) and (b) the tail of 
the bunch is at t = 0. 


A. Case of an ideal laser-shaping technique 

We base our approach on Ref. [21] where we developed 
a simple ID space charge model to investigate the expan¬ 
sion forces in various distributions; for polynomial 
distributions, we observed that the relatively small fields 
in the bunch-head will essentially preserve the local lon¬ 
gitudinal form. In the rear portion of the bunch however, 
there is an asymmetrical blowout which has proper sign 
to possibly lead to a linear-like tail which is required for 
the quadratic-ramp. 

We consider a laser intensity distribution of the form 
/(r, t) = T(t)R(r), where T{t) is the temporal profile 
and R{r) the transverse laser envelope. In our previous 
studies we used a radially uniform transverse distribu¬ 
tion and explored polynomial and exponential forms for 
T{t). In this section we report on the performance of the 
polynomial distribution given by 

T{t)=Tot^H{r-t), (17) 

where Tq is a normalization constant, a > 0 is the poly¬ 
nomial power, r is the ending time of the pulse, and H{t) 
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TABLE L Table comparing several different proposed drive bunch distributions as a function of bunch length and charge. 
Additionally, the maximum decelerating held (ET) is shown for each distribution. Here we consider k. — 1. 


distribution 

TZ{N) 

TliQ) 

FT 

doorstep [8] 

+ (1 - 7r/2 + 27rAf)2 

^2 + ni^-l) 

aX 

TT 

double triangle [9] 

yi + ( 27 riV - 1)2 

^2 + nC^-l) 

aX 

TT 

sin {q = 2) 

|y7r2(3- 16iV)2 +64 

yQ4-15n^ + ^ 

16aA 

Stt 

sin {q = 3) 

|y 7 r 2 (l -4iV)2 +4 

i^44-97r2 + 2^ 

6aX 

TT 

quadratic 

yi+ 7 r 2 ( 2 iV-l )2 


aX^ 

TT^ 


is the Heaviside function. 

The exponent a greatly influences the space-charge 
helds. Large values of a (e.g a > 5) lead to large space- 
charge forces which results in a uniformly hlled ellip¬ 
soidal distribution [22]. Alternatively, smaller values of 
a {a ^ 3) lead to a more uniform evolution of the bunch 
dynamics due to the increased uniformity of the held over 
the bunch. Additionally, the transverse spot size of the 
laser pulse on the photocathode also controls the longi¬ 
tudinal electric helds but also inhuences the transverse 
“thermal” emittances. It is also possible to reduce the 
electric helds and the associated blowout rate by using 
longer laser pulses; in this scenario, the resulting electron 
bunch will evolve at a slower rate but the resulting bunch 
distribution will have a smaller peak current compared to 
when starting with smaller values of r. A smaller current 
will impact the performances of the wakeheld accelerator 
(or require the implementation of a longitudinal compres¬ 
sion scheme). Finally, it would also be possible to use a 
longer, e.g. 2 + ^-cell, RF gun or another acceleration 
cavity in close proximity to the gun to preserve larger 
charge densities which could effectively alleviate the need 
for a bunch compressor to drive large accelerating helds 
in the subsequent wakeheld accelerator. 

Figure 5 shows simulated longitudinal phase space 
snapshots and corresponding currents at different axial 
locations downstream of the gun for a 1-nC bunch. For 
this simulation a 1-mm rms laser spot size on the pho¬ 
tocathode was used. The initial laser distribution was 
described by Eq. 17 with a = 2 and r = 15 ps. A ht 
of the current distribution at s = 50 cm from the photo¬ 
cathode is shown in Fig 5 and indicates that the hnal elec¬ 
tron bunch distribution is indeed accurately described by 
Eq. 13. 


B. Limitation of a practical laser shaping technique 

As a hrst step toward a realistic model for the achiev¬ 
able shaped we consider the photoemission process to 
be resulting from frequency tripling of a Aq = 800-nm 
amplifier infrared (IR) pulse impinging a fast-response 
time cathode (with typical work functions corresponding 
to ultraviolet photon energy ^ Ao/ 3 ). Such a setup is 




FIG. 5. Evolution of the electron-bunch current (a) and lon¬ 
gitudinal phase space (b) along the beamline at 20 (red), 60 
(green), and 100 cm (blue) from the photocathode surface 
and (c) comparison of the current profile numerically simu¬ 
lated at s = 50 cm (red trace) with a fit to equation Eq. 13 
(blue lines). The head of the bunch is at large values of z. 

commonly used in RF photoinjectors such as the one 
discussed in the previous sections. We further assume 
that the frequency up-conversion process does not affect 
the original laser’s temporal shape (e.g. the UV-pulse 
temporal shape is identical to the IR-pulse temporal 
shape). Under such an assumption, the formation of the 
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ideal temporal shape discussed in the previous Section 
is limited by the finite laser bandwidth and frequency 
response of the shaping process. 

We consider an incoming amplified IR pulse with inten¬ 
sity linir^t) = /o(r)sech^(t/r) downstream of the last- 
stage amplification, where r is the laser pulse duration. 
We model the IR pulse laser-shaping process via the 
convolution where 

dout{t) and R{t) represent the shaped-pulse intensity and 
response function of the shaping method respectively. 

Given the desired output shape and incoming laser 
pulse profile, the response function of the shaping pro¬ 
cess has to be set to satisfy [23] 

RH = (18) 

where the upper tilde represents the Fourier transfor¬ 
mation f{uj) = In practice Iin{^) is de¬ 

fined over a finite range of frequency uu = ujq ^ ^ where 
ujQ = ^ is the central laser frequency and 5uj = is 
the laser pulse bandwidth (^A is the wavelength span of 
the pulse spectrum). 

The typical shape considered in the previous section 
after laser shaping is shown in Fig. 6; the limited band¬ 
width has very little effect except for the well-known ring¬ 
ing effect at the sharp discontinuities [24]; see Fig. 6 (b) 
and (c). 

Another potential limitation to our shaping scheme 



FIG. 6. Comparison of nominal (“ideal”) shape with the 
shapes achieved when taking into account the photoemis¬ 
sion response time (“cathode”), the laser-pulse-shaping fi¬ 
nite bandwidth (“shaping”) and both effects (“cathode -h 
shaper”). The ideal laser temporal profile is described by 
Eq. 17 with 0 = 2 and r = 15 ps. Insets (b) and (c) are 
zooms of the areas t G [—15200, —13600] fs (peak location) 
and t G [—16000, —15020] fs (left edge of the profile) respec¬ 
tively. The head of the laser pulse is at t = 0. 

arises with a high-efhciency (semiconductor) photocath¬ 


ode. We consider as an example the case of Cs 2 Te photo¬ 
cathodes because of their wide use in high-current pho¬ 
toinjectors. The response-time limitation is investigated 
using the parameterized impulsional time response of 
Gs 2 Te described in Ref. [25] based on numerical simu¬ 
lations presented in Ref. [26]. The impulsional response 
is convolved with the distribution used in the previous 
section and the results are gathered in Fig. 6. Again this 
effect appears to be marginal. For the sake of complete¬ 
ness, the various profiles shown in Fig. 6 are tracked with 
ASTRA and the final current distributions at s = 50 cm 
are found to be indiscernibly close to the ideal shape 
considered in the previous Section; see Fig. 7. Such a 
result gives further confidence in the proposed shaping 
approach. 



FIG. 7. Comparison of the final electron-bunch current at 
s = 50 cm from the cathode surface for the four cases con¬ 
sidered in Fig. 6. The “cathode” and “shaper” respectively 
correspond to the inclusion of the cathode response time and 
shaper bandwidth limitation in the initial particle distribu¬ 
tion at s = 0 while the ideal case is given by by Eq. 17 with 
0 = 2 and r = 15 ps. The head of the bunch corresponds to 
z > 0. 


IV. FORMATION OF HIGH-ENERGY 
TAILORED BUNCHES FOR A DWFA LINAC 

We finally investigate the combination of the tailored 
current-profile generation scheme with subsequent accel¬ 
eration in a linac located downstream of the RF gun. 
Such a configuration could be useful to form tailored rela¬ 
tivistic electron bunches for direct injection in wakefield- 
acceleration structures. For this example, we consider 
a high-repetition drive bunch with parameters consis¬ 
tent with a recently proposed beam-driven accelerator 
for a short-wavelength free-electron laser (FEE) [14]. We 
adopt a different approach than Ref [14] and instead 
choose a 1.3-GHz superconducting RF (SGRF) linac (LO 
and LI) composed of TESLA cavities [29] coupled to a 
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quarter-wave 200-MHz SCRF gun [30, 31] originally de¬ 
signed for the WiFEL project [32]; see diagram in Fig. 8. 
The accelerator also includes a 3.9-GHz accelerating cav¬ 
ity (L39) section to remove nonlinearities in the longitu¬ 
dinal phase space [33, 34]. For this study we explored the 
use of polynomial laser profile described by Eq. 17 and 
let a and r as free parameters. 



; ID longitudinal-dynamics 


ASTRA program 

FIG. 8. Block diagram of the accelerator configuration ex¬ 
plored for the formation of high-energy ramped bunches. The 
legend is as follows: “QW” stands for quarter-wave, “LO” 
and “LI” are standard 1.3-GHz cryomodule equipped with 8 
TESLA-type SCRF cavities, “L39” is a cryomodule consist¬ 
ing of four 3.9-GHz cavities, and ”BC” is a magnetic bunch 
compressor. 

The laser-profile parameters and accelerator settings 
were optimized using a genetic optimizer [35] to result 
in a final distribution with current profile consistent to 
achieve a high transformer ratio. The optimized acceler¬ 
ator settings are summarized in Tab. II. In our optimiza¬ 
tion, we chose the wakefield structure to be a dielectric- 
lined waveguide with parameters tabulated in Tab. Ill 
and we introduce a longitudinal scaling factor r] as free 
parameter such that the axial coordinate is scaled follow¬ 
ing z ^ z' = rjz. The optimization converged to a value 
T] = 0.16. The obtained wakefield and scaled shape are 
shown in Fig. 9 (a). For the wakefield calculations we 
followed the formalism detailed in Ref. [39] and use the 
first four modes in the wakepotential used for the beam 
dynamics simulations. 

Given the devised configuration, a one-dimensional 
model of the longitudinal beam dynamics was employed 
to asses the viability of the required compression and 
especially explore the possible impact of nonlinearities 
in the longitudinal phase space on the achieved cur¬ 
rent profile. We considered the current could be lon¬ 
gitudinally compressed using a conventional magnetic 
bunch compressor (BG) with longitudinal linear and sec¬ 
ond order dispersions R^e and T 566 = — §^56 [36]. In 
our simulations the longitudinal dispersion was taken 
to R 56 = —20 cm following similar designs [37]. The 
phase of LO and phase and amplitude of L39 were em¬ 
pirically optimized and the resulting longitudinal phase 
space {zq^ ^ 0 ) was tracked through the BG via the trans¬ 
formation zq ^ z = zq -\-T^qqSq . An optimum set 

of phases and amplitudes was found and listed in Tab. II 
and the sequence of the longitudinal phase spaces along 
the injector appear in Fig. 10. The final wakefield excited 
in the structure with parameters listed in Table III is dis¬ 


TABLE 11. Optimized settings for the accelerator param¬ 
eters needed to produce and accelerate a ramp bunch to 
~ 200 MeV. The parameter a and r are defined in Eq. 17. 


parameter 

value 

units 

laser rms spot size cr^ 

2.5 

mm 

laser ramp a parameter 

19.86 

- 

laser ramp duration r 

96.8 

ps 

bunch charge Q 

5 

nC 

peak E-field on cathode 

40 

MV/m 

laser injection phase 

71.0 

deg (200 MHz) 

gun output beam momentum 

5.15 

MeV/c 

ace. voltage LO 

165 

MV/m 

off-crest phase LO 

-12.35 

deg (1.3 GHz) 

ace. voltage L39 

24.1 

MV 

off-crest phase L39 

-192.35 

deg (3.9 GHz) 

beam momentum after L39 

- 143 

MeV/c 

final beam momentum after LI 

- 350 

MeV/c 


3.5, -^^^^^^-^80 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 


z (mm) 

FIG. 9. Final current distribution (green shaded area) and 
associated wakefield (blue traces) for the “ideal” (a) and “re¬ 
alistic” (b) cases of compression discussed in the text. The 
head of the bunch corresponds to z = 0 


played in Fig. 9 (b) — the achieved field and transformer 
ratio values are summarized in Table III. We remark that 
the inclusion of a refined model of longitudinal dynamics 
leads to the apparition of features [e.g. a small current 
spike in the bunch tail; see Fig. 9 (b) or 10 (d)] that were 
absent in the optimization process implementing a sim¬ 
ple scaling of the longitudinal coordinates; see Fig. 9 (a). 
The origin of the small current spike can be traced back 
to the nonlinear correlation imposed by space charge in 
the early stages of the bunch-transport process (i.e. in 
the drift space upstream of LO); see Fig. 10 (a). Never¬ 
theless the achieved peak field and transformer ratio as 
the bunch passes through the DLW are very close (within 
10%) to the ones obtained with the scaled distribution. 
These results indicate that our proposed injector concept 
appears to produce the required current profile. Further 

















- 0.8 

- 1.6 


0.0 0.5 1.0 1.5 



-2.5 0.0 2.5 5.0 
2: (mm) 



-2.5 0.0 2.5 5.0 
2: (mm) 



-0.60.0 0.6 1.2 
2: (mm) 


FIG. 10. Snapshots of the longitudinal phase spaces and as¬ 
sociate current profiles (red traces) upstream of LO (a) and 
downstream of LO (b), L39 (c) and BC (d). Simulations up 
to L39 are carried with ASTRA whereas a one-dimensional 
longitudinal-dynamics model is used for BC2. The head of 
the bunch corresponds to 2 > 0. 


studies, including a transverse beam dynamics optimiza¬ 
tion and the inclusion of collective effects such as coher¬ 
ent synchrotron radiation and space charge downstream 
of L39 and throughout the bunch compressor, will be 
needed to formulate a detailed design of the injector. We 
nevertheless stress that the simple model presented above 
confirms a plausible longitudinal-beam-dynamics capable 
of preserving the formed current profiles after accelera¬ 
tion and compression. The final energies and peak cur¬ 
rents are all within the parameters suggested in Ref. [14]. 


TABLE III. Dielectric-line waveguide (DLW) parameters and 
resulting wakefield values using the current profile shown in 
Fig 9. The “ideal-” and “realistic-compression” entries re¬ 
spectively correspond to the cases when the final current pro¬ 
file is obtained via a simple longitudinal-axis scaling or via 
particle tracking. 


parameter, symbol 

value 

units 

DLW inner radius, 

750 

fim 

DLW outer radius, Vo 

795 

/j.m 

DLW relative permittivity, Cr 

5.7 

- 

DLW fundamental mode, fi 

369.3 

GHz 


ideal compression: 

Peak decelerating field, \E-\ 14.01 MV/m 

Peak accelerating field, |E+| 75.55 MV/m 

transformer ratio, 77. 5.39 

realistic compression: 

Peak decelerating field, \E-\ 12.84 MV/m 

Peak accelerating field, |E+| 63.87 MV/m 

transformer ratio, 77. 4.95 


We finally note that the generated current profiles are ca¬ 
pable of supporting electric fields and transformer-ratios 
in a DLW structure with performances that strike a bal¬ 
ance between the two cases listed as “case 1” and “case 
2” in Table 1 of Ref. [14]; see Tab. III. A simple estimate 
indicates that our drive bunch would require a DWFA 
linac of ^ 30 m in order to accelerate an incoming 350- 
MeV witness bunch to a final energy of ^ 2 GeV. 


V. SUMMARY 

In conclusion, we have presented a set of smooth 
current profiles for beam-driven acceleration which 
display comparable performances with more complex 
discontinuous shapes discussed in previous work. We 
find that all proposed current profiles which lead to 
uniform decelerating fields “live” on the same perfor¬ 
mance curve and that a given profile can be scaled 
to a particular accelerating field or transformer ratio. 
We also presented a simple laser-shaping technique 
combined with a photoinjector to generate our proposed 
quadratic current profile. We finally illustrated the 
possible use of this technique to form an electron bunch 
with a tailored current profile. The distribution obtained 
from these start-to-end simulations were shown to result 
in a transformer ratio ^ 5 and peak accelerating field of 
E_^ ^ 60 MV/m in a dielectric-lined waveguide consis¬ 
tent with the proposal of Ref. [14]. The method offers 
greater simplicity over other proposed techniques, e.g., 
based on complex phase-space manipulations [11, 38]. 
Finally, we point out that the proposed method could 
provide bunch shapes consistent with those required to 
mitigate energy-spread and transverse emittance dilu¬ 
tions due coherent-synchrotron-radiation in magnetic 
bunch compressors [40]. 
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Appendix A: Maximum of Acos{kz) + B sm{kz) 

The accelerating field behind the bunch often assumes 
the functional form 

F{z) = Acos{kz) + Bsm{kz). (Al) 

The procedure to evaluate the transformer ratio entails 
to determining the maximum value of F{z). Such a value 
if found by solving for 

= k[—A sm{kz) + B cos{kz)] = 0, (A 2 ) 

with solution 2^5 given by 

tdin{kzs) = — = T. (A3) 

Squaring the previous equation, it is straightforward to 
show that 

1 

sin2(A:2:^) = cos^{kz^) = (A4) 

Expressing the value oi F{zs) using the previous equa¬ 
tion in Al leads to the maximum value of F{z) 


In this Appendix we summarize and rewrite in nota¬ 
tions consistent with our Section II the equations de¬ 
scribing the linear ramp [ 8 ] and double-triangle [9] cur¬ 
rent profiles. These equations are the ones used in Sec¬ 
tion IIC. 

The “doorstep” current profile considered in Ref. [ 8 ] is 
written as 




a if 0 < ^ < ^, 

< + 1) \ii<z<Z, 

0 elsewhere. 


(Bl) 


The “double-triangle” suggested in Ref. [9] is given in 
our notations as 


S'( 2 :) 


akz 

< a{kz — 1 ) 

0 


if 0 < 2 ; < (^, 
if^<z<Z, 
elsewhere. 


(B2) 


F = F{zs) = A^l+T^. (A5) 

The latter equation is used at several instances through¬ 
out Section II. 


For both cases, ^ = A/4 leads to the fiat decelerating 
fields over the tail of distribution and leads to the FT 
and IZ tabulated and illustrated in Tab. I and Fig. 3 
respectively. 
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